MATRIX PERMANENT AND QUANTUM ENTANGLEMENT 
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Abstract. We point out that a geometric measure of quantum entanglement is related to the 
matrix permanent when restricted to permutation invariant states. This connection allows us 
to interpret the permanent as an angle between vectors. By employing a recently introduced 
permanent inequality [Carlen, Loss and Lieb, Meth. and Appl. of Analysis, 13 (2006), no. 
1, 1—17], we can write combinatorial formulas for quantifying the entanglement of permutation 
invariant basis states. When applying the geometric measure to permutation invariant states 
with nonnegative coefficients, we show that the overlap with a product state is maximized by a 
tensor product of the same single-party state. This extends some observations in [Hayashi et al., 
Phys. Rev. A 77, 012104 (2008)]. 

1. Introduction 

In the Editor's statement forewording the 1982 monograph Permanents by Mine [28], Gian-Carlo 

Rota wrote the following words: 

"A permanent is an improbable construction to which we might have given little 
chance of survival fifty years ago. Yet numerous appearances it has made in physics 
and in probability betoken the mystifying usefulness of the concept, which has a 
way of recurring in the most disparate circumstances." 
The present paper highlights a connection between the permanent and entanglement of certain 

quantum states, therefore putting in evidence a further appearance of the permanent in physics. 

Some background is useful for delineating the context. The permanent and the determinant of an 

n x n matrix A (with entries in a commutative ring) are respectively defined as 

n n 

perm(A)= £ ]jA lMl) and det (A) = £ fi^i), 

■rr£S n 8=1 7reS„ i=l 

where S n denotes the full symmetric group on a set of n symbols. In light of the definitions, it 
is reasonable that there are cases, as it was originarily observed in 1859 by Cayley [7j, in which 
permanents can be computed by means of determinants (see also, e.g., Kasteleyn [22], Godsil and 
Gutman [T2], and Frieze and Jerrum [H]). Geometrically, the determinant is the volume of the 
parallelepiped defined by the lines of the matrix; algebraically, it is the product of all eigenvalues 
including their multiplicities. It is curious that despite its striking similarity with the determinant, 
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the permanent does not have any known geometric or algebraic interpretation. Moreover, while 
standard Gaussian elimination provides an efficient technique for computing the determinant, the 
exact computation of the permanent remains a notoriously difficult problem. 

The best known algorithm for an nxn matrix, due to Ryser in 1963 [33], needs Q(n2 n ) operations. 
By a seminal result of Valiant [38] , computing the permanent is indeed "The #P-hard problem" . 
Thus computing the permanent on worst case inputs cannot be done in polynomial time unless 
P# p = P and in particular P = NP. It follows that algorithms for the permanent acquired a special 
position in computational complexity to the point of becoming a fertile ground for many approaches, 
including iterative balancing (Linial et al. |26p. elementary recursive algorithms (Rasmussen }32] ) 
and, most relevantly, Markov chain Monte Carlo methods (Broder 0], Jerrum and Sinclair [21] )• 
These efforts contributed to a deeper understanding of the permanent and produced important 
results as fully polynomial randomized approximation schemes for several types of matrices (see 
Jerrum, Sinclair and Vigoda [25], Barvinok [2], and the references contained therein). 

From the mathematical perspective, there are two main lines of research centered on the perma- 
nent: the study of permanents and probability (see, e.g., Friedland et al. [10]); min/max questions 
concerning a number of inequalities, peaking with the proofs by Egorychev [8J and Falikman [S] 
(see also [3D]) of the 1926 van der Waerden's conjecture about permanents of doubly stochastic 
matrices. 

When focusing our attention on physics, we may divide into at least four groups the known 
applications of the permanent: the many aspects of the dimer problem, some uses involving the 
hafnian, Monte Carlo generators, and, more pertinent to our discussion, linear optical networks for 
quantum information processing. Intuitively, the permanent tends to be related to bosons while 
the determinant to fermions. A list of applications follows. 

The problem of computing the permanent of a (0, l)-matrix is the same as the problem of counting 
the number of perfect matchings in a bipartite graph. This translates into the dimer problem 
[22l [T9] . The problem is traditionally related to models for adsorption of diatomic molecules on 
crystal surfaces, mixtures of molecules of different sizes and the cell-cluster theory of the liquid state 
(see, e.g., Welsh [42]). New applications occur in the study of configurations of melting crystals 
(Okounkov et al. [27]), BPS black holes (Heckman and Vafa [18]), and quiver gauge theories 
(Hanany and Kennaway [14] ) . 

The hafnian of a matrix A, denoted by hf(^4), is a polynomial generalizing the permanent as 
the pfaffian generalizes the determinant. This notion was originarily introduced by Caianiello to 
express the perturbation expansions of boson field theories [5] . 

Bialas and Krzywicki [3] (see also Wosiek [43]) introduced a procedure to include Bose-Einstein 
correlations in Monte Carlo event generators. The procedure makes use of the generalized Wigner 
functions and it requires to compute the permanent of a correlation matrix depending on particle 
momenta and other parameters. 

Sheele et al. |35j have proved that matrix elements of unitarily transformed photonic multi-mode 
states can be written as permanents associated with the symmetric tensor power of the beam splitter 
matrix (see also Kok et al. 23J). The result implies that computing matrix elements in the Fock 
basis is not an easy task. Permanents count the ways of redistributing n single photons through 
an SU(n) network to yield exactly n single photons at the outputs and then allow to compute 
probability amplitudes [3H [3S] ; optimal networks are obtained by maximizing the permanent under 
given constraints. 

The present paper is concerned with quantum entanglement. It is now established that entangle- 
ment is an important physical quantity whose presence appears to be necessary in many applications 
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of quantum information processing [31] . Given a generic quantum state, detecting and measuring 
its entanglement is a challenge from both the mathematical and the experimental point of view. 
It is computationally hard [13] and furthermore there is no general agreement on how to quantify 
entanglement. A number of different entanglement measures have been therefore introduced using 
a variety of approaches (see the Horodeckis' review [20] and its references). 

Here we consider a geometric measure. This quantifies the angle between two states, subject to 
an optimization problem: the closest state with zero entanglement and the state under analysis. 
It turns out that, when we restrict the analysis to permutation invariant states, such a measure 
can be described in terms of the permanent of a certain matrix. A recent permanent inequality of 
Carlen, Loss and Lieb [5] can be employed to bypass the optimization problem. Hence, maximizing 
the permanent of a matrix with certain constraints, reduces to a simpler problem. As we have 
mentioned earlier, the permanent does not have any known geometric interpretation. The present 
work contributes in this direction, by interpreting the permanent as an angle between vectors (or 
a cosine, to be more precise). Similar results concerning entanglement of permutation invariant 
states are obtained independently by Hayashi et al. [17] without the use of permanents. 

The rest of the paper is structured as follows. In Section 2, we define the geometric measure of 
entanglement. In Section 3, we state and prove our results. For the sake of clarity, we shall also 
include a few examples 



2. Geometric measure of entanglement 

The geometric measure of entanglement used here was firstly introduced by Shimony [37] in the 
setting of bipartite pure states. It was generalized to the multipartite states by Barnum and Linden 
P] , and further extended by Wei and Goldbart [41] . The intuition beyond the measure consists of 
thinking about entanglement as an angle between two states: namely, the state in analysis and a 
product state, i.e., a state with zero entanglement. Crucially, the product state is chosen over all 
possible product states so that it minimizes the angle with the state in analysis. 

Let 



=(l) e (2)... e (n)\ 

"PI C P2 C p™ / 



(2-i) m= E w-J^M 

Pl'-'Pn 

be a generic multipartite pure state in a Hilbert space H. = ® k of dimensionality dim(7i) = 

n^=i h- Each set {\e p 1 }) : p t — 1, 2, hi} is a local basis for the i-th subspace C^. A pure state 
in H. is said to be a product state if it can be written in the form 

n n 

i0) = (g)i^ i) ) = (g)E( c S ) i^))' 

i— 1 i— 1 Pi 

where \<j>^) €E C^ 4 is some pure state; \ip) is said to be entangled, otherwise. Given a state \ip) as 
in Eq. (|2.1[) . let us define 

A m axO0) := max|(</>|V>)|, 
\4>> 

where the maximization is performed over all product states \<p) This formula tells us how well 

the possibly entangled state can be approximated by a product state. The formula provides a 
method that satisfies the desiderata for a well-defined measure of entanglement [20] . Such a method 
is usually called geometric measure. The terminology is justified since A max (tp) is an angle between 
two vectors. Hence, notice that the amount of entanglement increases while A max ( , 0) decreases and 
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therefore the quantity of entanglement depends essentially on A max (tp). Given a state concrete 
geometric measures are 

E sin 2 (ip) := 1 - A^^) and E los (ip) := -2 log 2 A max (», 

introduced in [41] and [40], respectively. The relation of the geometric measure to other measures 
has been studied in [15l [15l EHJ • In the next section, we will show how the geometric measure can 
be related to the permanent. 

3. Permanent and entanglement 
A permutation invariant basis state is a pure state of the form 



\S(n,k)) 



MV^V^,..,^..^)), where C£ := n!/ JJ hi 
- _ i=i 



ki k 2 k d 

As an example, we consider the permutation invariant basis state 
V6 



\S(4, (2, 2))) = ^ E 1^(1,1,2,2)) = -L (|H22) + |1212) + |2112) + |1221) + |2121) + |2211)) . 



V6 



Here n = 4 and d = 2. Theorem [T] is our main result concerning these class of states: 
Theorem 1. Let \S(n, k)) be a permutation invariant basis state. Then 



( s <^Wtftf n 

V lli=l "-*• i=l;ki^0 



Proof. We prove the statement by comparing the possibly entangled state \S(n, k)) to the general 
product states 

n / d \ d 

i$ = <8) with Ek<i 2 = l 

j=i ) i=i 

According to the definition of geometric measure, the first step is to evaluate the overlap 



which gives 
(3-1) 0g 



4>,:= (S(n;k)\<f>), 



C'l JC? 

2_j ^tC 1 ). 1 ■ ■ ■ a 7r«(fei),l a 7r i (fci+l),2 ■ ■ ■ a m(n),d = — P er (^fc) ; 

7T 4 eS„ U - 



where is annxn matrix defined as follows: 
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The matrix has fcj identical columns vi = (ai,i, ct2,i, ■ ■ ■ > a n ,i) T with i = 1, . . . , d and X)j=i ^ = 
n. The next step consists of maximizing the absolute value of the overlap, i.e., \4>%\, over the set of 
all otj/s such that Yli=i \ a j.i\ 2 = 1- O n the basis of Eq. ()3.1|) . we can write 



A max ( fc) = max — per(A £ ) 



To deal with this equation, we will use a recent result by Carlen, Loss and Lieb [6] . For any matrix 
F defined as 

F ■= [/l, /2, ■ ■ ■ , fn\) 

where ft, fa, . . . , f n are arbitrary column vectors of dimension n, they have shown that 



(3.2) 



i=l 



where | |/i 1 12 denotes the L2-norm. The r.h.s. of the inequality can also be regarded as the permanent 
of a matrix whose i-th column contains only identical entries \fi\/^/n . For example, 



\\fi\J 



f\f 2 \\ 
1/1 

\\f2\J 



f\fn\\ 

l/n| 

V|/n|/ 



By applying the inequality in Eq. (|3.2[) . we obtain 



where 

(3.4) ai 
with the property that X^j=i a i = ^ an d 



\ 



1 " 

n — ^ 



3=1 
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(a 2 \ 



a 2 



( a2 \ 

Ot-2 



\aij \a 2 / \ a 2/ 
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By Eq. l|3~5]) . we have 

n 

A max I S(n, k) ) — max </>,-< max wCsTTa^ = vaax.{<bq\S(n: k)), 
V V U) k a,eR+ V k 11 ' <t> a 
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where 
(3.5) 



10s) 



3=1 \Z=1 



The interpretation is that there is a product state constructed from \<p) by appropriately averaging 
the coefficients as in Eq. p. 41) and that the derived product state has a larger overlap. The resulting 
\4> s ) is a tensor product of n copies of the same state. The last step is a simple maximization 
procedure. We need to maximize the function 



/ (Xl,X 2 , :.,X d ) 



with nonnegative domain, under the constraint x i = 1* This gives 



A M ,(5(n,g))= > /q J] (~ 



which verifies the statement. 



□ 



Note that in order to find A max (|5(n, fe))), it is sufficient to use the state \4>) to be the product 
of n identical copies of an arbitrary single-party state \a), i.e., 

(3.6) \4>)=® 7 UW). 

When the number of levels d is equal to the number of parties n and fcj = 1, for every i = 1,2, n, 
we have 



The form of the matrix is generic. The only constraint is that each column has unit norm. In 
this case, 

















«2,1 


«2,2 




OL2,d 






\ a d,lJ 


\ a d,2/ 









> d u\ 



A max (5(d, (1,1, ...,!))) 



For the case of qubits, namely when d = 2, we can prove the corresponding result without using 
the Carlen-Lieb-Loss inequality, but the older Schwarz and McClaurin inequalities. In this case, 
the permutation invariant states have the form 

\S(n,k)) = ^J2^ eS Jn i (0^,h : A)), where C? := j^g^. 

k n — k 

The theorem below states the connected result: 

Theorem 2. Let \S(n,k)) be a permutation invariant basis state for qubits. Then 

n—k 

' n — 



A max («, fe) 



fc! (n — fe)! \n 



PERMANENT AND ENTANGLEMENT 



7 



Proof. We want to find the maximal overlap between \S(n, k)) with product states 
As the coefficients in \S(n, k)) are nonnegative, we can set /3j = 0. We then evaluate 



/ (J 71 I ''' 



<t>k = {S{n,k)\4>) = — p ^ IIv^W II V 1 -^i(fe+i) 



ra! 

7T i^zS n \l = l l = k + l 



Using the Cauchy-Schwarz inequality, we have 



(n!) 

By the Maclaurin inequality 



i^i 2 ^7^( e e n i-^oH-D 



1 fe /l " 

^ E IKw< -E 

7ri6S„i=l \ i=l 



V 



we arrive at 

l^fel 2 < C fc n (g) fe (l - q) n - k = Ct cos 2fc flsin 2 ^ 0, 
for cos 2 8 = q. This means that 

141 < V^(f) fc/2 (1 - q) (n ~ k)/2 = v/Cfcos^sinf"-^ 0. 
Maximizing the expression on the r.h.s. over the angle 9, we obtain 
A max (n, k) = max v /C^cos fc 6 sin ( "" fc) 8 



k \ 2 / n — k 



y fc! (n — fc)! \n j 

This concludes the proof. □ 

For the case of three qubits, some examples of permutation invariant states are the following 
ones: 

|GHZ> = (|000) + |111))/V2, |W) = (|001) + |010) + |100))/V3, |W) = (|H0) + |101) + |Qll))/\/3. 

The states |GHZ) and |W) have extremal properties and have particularly important roles in quan- 
tum mechanics [20]: the Greenberger-Horne-Zeilinger state, |GHZ), was used to test Bell's inequal- 
ities; the W state, |W), exhibits genuine three-party entanglement, in a different way from |GHZ). 
By applying Theorem^ we have 

A raax (GHZ) = l/y/2 and A max (W) = A max (W) = 2/3. 

Let us now consider examples for three-party 4-level systems, namely n = 3 and d = 4. The 
chosen vectors are a — (2, 0, 0, 1) and b — (1, 1, 1, 0) and their corresponding states are 

\a) = ^(|H4) + [141) + |411>), and \b) = ^(|123) + |132) + |213) + |231) + |312) + |321>). 

With the use of Theorem [1] we have 

A max (a) = A max (W) = 2/3 and A max (6) = y/2/3. 

Note that the states \a) and |W) have the same structure. 

We conclude by asking the following question and then providing a partial answer: 
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Problem 1. Is it true that in order to obtain the maximal overlap of any permutation invariant 
state, we can assume the product state to be a tensor product of the same single-party state ? 

Hayashi et al. [16] have attempted to answer this question. However, it still remains an open 
problem. Here we prove that the question has an affirmative answer for certain classes of permuta- 
tion invariant states, beyond the basis states discussed above. 

The first class we consider is the case where the coefficients c^'s are nonnegative. To obtain 
the maximal overlap for the corresponding state l^), we can as well set the coefficients in the 
unentangled state \cf)) to be nonnegative, as the goal is to maximize the overlap between |^>) and 
\(f>). Thus we have 

m - E = E ^-^p^) ^ E c s ^r^< A k) - E c * ■ *5 4 = <^>, 

where we have used again the inequality proved in [BJ. This means that for nonnegative cg's, in 
order to maximize the overlap, we can use the product state \<p s ) in Eq. (|3.5p . consisting of a direct 
product of identical single-party states. An example of this is given by the states 

| WW (s)) = </a\ W> + A /T _ s| W) , where s e [0, 1] . 

The way to use the product state in the form of Eq. (|3.6|) is justified. This was the ansatz used to 
calculate the entanglement for this family of states in [IT] . In particular, the product state can be 
written as 

\<f> 8 (9)) = (cos6>|0) +sin6»|l)) 83 

and then, by maximizing the inner product (<p s {9) |WW (s)) over 9, it is straighforward to obtain 
the expression 

A max (WW (s)) = - (V^cos 9(s) + VT~s~ sin 9 (s)) sin 29(s), 

where 9 (s) is the solution of the equation 

VT - 7tan 3 6» + 2 % /itan 2 6»-2\/T _7 7tan6'- Vs = 0, where tan0 e [l/V%, V%\ . 

We can also approach a more general class of states. When the coefficient c^' s are arbitrary 
but nonnegative, the above consideration of using states as in Eq. (|3 . 6[) holds; this gives us the 
corresponding state 

\^) = ^ % \S(n,k)). 

k 

Now, we perform a basis change on \ip), i.e., 

= {u ® u ® ... ® u)\i>) = J2 h\ s ^ *)>. 

k 

where U is any unitary transformation in U(d). Specifically, the transformation U acts on a single 
party, and b^s are the resulting coefficients for \ip r ) expanded in the basis of \S(n, k)}. The resulting 
coefficients 6^'s are in general complex. Since we have shown that to calculate the entanglement for 
\ip) we can assume the product state to be as in Eq. (|3.6[) and \tp ) is simply given by a local change 
of basis, to calculate the entanglement for l^'), we can take a fiducial state of the same form. To 
illustrate this fact, we consider a generic element of SU(2): 

U=( U , " ], where \u\ 2 + \v\ 2 = 1. 
V — v u J 
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The effect of U on the parties of |W) is given by 




and similarly for |W), 1 000) and | 111), our basis states in the symmetric subspace. The correspond- 
ing coefficients are in general complex. When U is diagonal, the coeffients are transformed as 

tfjt = c^e lk ' e , where 9 is an arbitrary real d-component vector characterizing the matrix U. 

It takes 8 real parameters to describe the generic permutation invariant states for three qubits. 
If we start with 4 arbitrary nonnegative coefficients and supplement with arbitrary U(2) transfor- 
mations (which have 4 real parameters) , we have then 8 real parameters in total. The number boils 
down to 6 in both cases, if we take into account normalization and global phase. This counting 
suggests that the statement in Problem [1] may well be true for three qubits. 
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